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We obtain a general result for the Lamb shift of excited states of multi-level atoms in inhomoge- 

neous electromagnetic structures and apply it to study atomic hydrogen in inverse-opal photonic 

ff^ , crystals. We find that the photonic-crystal environment can lead to very large values of the Lamb 
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shift, as compared to the case of vacuum. We also predict that the position-dependent Lamb shift 

^: 

^^ ■ should extend from a single level to a mini-band for an assemble of atoms with random distribution 

^: 

^^ . in space, similar to the velocity-dependent Doppler effect in atomic/molecular gases. 



PACS numbers: 42.70.Qs, 32.80.-t, 42.50.Ct. 



Since the pioneering experiment performed by Lamb and Retherford |l| in 1947 and the 
subsequent theoretical analysis developed by Bethe [2], the study of the Lamb shift plays an 
unique role in quantum electrodynamics (QED) because it provides an excellent test of the 
QED theory by comparing its predictions with experimental observations J3l,l4|. Recently, 
many efforts have been devoted to the study of various physical effects associated with the 
Lamb shift lilalll- 

Photonic crystals (PCs) are a new type of optical material with a periodic dielectric 
structure J8[. They can pronouncedly modify the photonic density of state (DOS) and 
local DOS lead ing to novel quantum-optics phenomena J2| such as inhibition [ll| and co- 
herent control 111 of spontaneous emission, enhanced quantum interference effects Il2|, 
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non-Markovian effects [l3|,|l^, wide lifetime distribution |l5|, non-classic decay [l6[, and 
slop discontinuities in the power spectra 'l7|, etc. 

Strong suppression or enhancement of light emission by the PC environment is expected 
to modify the Lamb shift. However, very different predictions for the Lamb shift can be found 
in literature. The isotropic dispersion model [l^ predicts an anomalous Lamb shift and level 
splitting for multi-level atoms. For two-level atoms, the anisotropic model [isf suggests that 
;he Lamb shift should be much smaller than that in vacuum, while the pseudogap model 
201 1 predicts a change of the Lamb shift of the order of 15% compared to its vacuum value. 
At last, a direct extension of the Lamb shift formulism for multi-level atoms in vacuum to 
the ca. of PC. .ugge.. that the La.b .Mft chtle. ne.H„b,y fro. .t. vacu.. value Q. 

Motivated by previous controversial results, in this Letter we employ the Green's function 
formalism of the evolution operator to obtain a general result for the Lamb shift in PCs. We 
reveal that in an inhomogeneous electromagnetic environment the dominant contribution to 
the Lamb shift comes from emission of real photon, while the contribution from emission 
and reabsorption of virtual photon is negligible, in vast contrast with the case of free space 
where the virtual photon processes play a key role. The properties of the Lamb shift near the 
band gap are calculated numerically for an inverse opal PC. We find that the PC structure 
can lead to a giant Lamb shift, and the Lamb shift is sensitive to both the position of an 
atom in PCs and the transition frequency of the related excited level. 

We study the Lamb shift in PCs in the framework of nonrelativistic quantum theory. For 
an multi-level atom located at the position r in a perfect 3D PC without defects, Hamiltonian 
of the system can be presented in the form H = Hq + H\^t + -f^ct, where the term H^ stands 



for noninteracting Hamiltonian and the term Hi^i describes interaction between an atom 
and photons, so that 
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Hq + Hint = h K + ^ 5Z ^nka;^k'^"k H P ■ A(r) (1) 
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with A(r) = ^^j^(^/2eoi^nk)^^^[E„(k, r)a„k + H. C] being the quantized vector potential, 
the second-order term of the vector potential in Eq. (^ has been neglected, and H^t = 
5rn/m)p'^/2rn is a mass-renormalization counter-term for an electron of observable mass m 
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22| . The electromagnetic (EM) eigenmodes {cj^k; E„k(r)} in PCs can be found by the 
plane- wave expansion method j23| . 

We assume that an atom is excited initially, and it stays at the /-th energy level without 
a photon in the EM field, and denote | / >=| /, > and | F;^^ >=| j, l„,k > (i.e. the atom is 
at the level j and the EM field has a photon in the state nk) as the initial and final states of 
the system, respectively. The state vector of the system evolves according to the equation, 
I ^(t) >= U{t) I / >= Ci{t) I / > +Ej,„kCnk(^) I Kk >> with the initial conditions 
Cj(0) = 1 and C„k(0) = 0, where U{t) is the evolution operator. Applying the Green's 
function technique to the evolution operator, we obtain the Fourier transform Ci{u!) of Ci{t) 
in the form [24], 

Q{u) = ±-[Gr^{u)-GUu)] , (2) 

with G^{uj) = lim^_^o+ < I I G{z = cu ±ir]) \ I >, where G{z) is defined by the operator 
identity G(z)(z — H/h) = 1. Projecting this operator identity onto the one-photon Hilbert 
space J25| and noting that the nonvanishing matrix elements of Hint are < Fl■^^ \ Hint \ I >, 
we obtain the following analytic expression 

Oti^) = lim ^ : ; — ; r, (3) 



where F = Y,- aijg{r, u -Uj), A = '^j{aij/2TT){uj - ujj)l3{r, oo - Uj), and 

^(r, cu) = ^ V / dk\ E„k(r) 1^ 6{uj - a;„k), (4) 

/5(r, u - u,) = V r' 9ir,uj') 

Jq {u - iUj - uj')uj' 



Here Vpc is the PC volume, uj^ei = mc^ /h is the relativistic limit of the photon energy 
aij = e^lpijl"^ /Sirm^eQhc^ is the relative line width of the atomic radiation from the /-state 



to j-state in vacuum, and V stands for the principal value of the integral. In Eqs. (@]) and 
(0), we have considered a random orientation of p;j and include the mass-renormalization 
contribution, respectively |l8|, |l9|, |20|, |2l|, |2^. The function g{v,uo) is the local spectral 



response function (LSRF) proportional to the photon local DOS. 

Equations (0) and show that the radiative correction to the bound level / is determined 
by the expression 
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{uj-uJi) = ^^{uj-ujj)[3{r,uj-Uj). (6) 
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In the two dispersion models, |E„kfr)| = l/lpc(V /m ), then Eq. (p) just gives the results 
described by Eq. (6a) of Ref. [lS| provided we take / = 1. For a two-level atom with j = 0, 1, 
we note that an = due to pn = (j>ij = i{uji — ujj)mrij), and Eq. (jHl) can be simplified 
to Eq. (4.9) of Ref. 20| by setting / = 1 and luq as zero point of energy. In vacuum, 
g{r, uj') = uj', and by setting u = ui in the right-hand side of Eq. (0), we obtain 

g2 

^? = ^22^3 Z2 -^.dP^il'/5(r, -uj^i), (7) 
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where Uji = ujj — uji and /3(r, —ojji) = — lia{ujrei/ujji + 1) ~ — ln{ujrei/ujji). Because /5(r, —ujji) 
is a slowly varying function of o^;, it is reasonable to make the approximation, Uji ^ u — ui, 
for (3{r,—iUji) (see also Ref. |22|), with u :^ lui being a weighted average of {ujj}. This 
approach implies that the dominant contributions to the Lamb shift come from the emission 
and reabsorption of virtual photons (corresponding to the transition processes from the 
I level to higher levels), rather than emission of real photon (corresponding to transition 
processes from the / level to lower levels). Noticing that '^j^ji\pij\'^ = ^e^|V^z(0)P/2£:o, 
where ipi{0) is the wave function value at the center of an atom in the state | / >, we finally 
obtain a standard nonrelativistic result, 

^' - 127r2m%gc3 ^^ C^ - CO,) ' ^^^ 

Thus, Eq. (0) gives a general result for a nonrelativistic radiative correction to a bound level 
of a multi-level atom in an inhomogeneous EM system. 

We solve Eq. (0 numerically for an actual PC structure. For calculating the function 
g{r, to'), we employ an efficient numerical method recently developed in Ref. 26|. For calcu- 
lating P{r,io — iOj), we make a reasonable approximation following the Refs. J21| and J27|: 
the dispersion function g{r,uj') of a PC vanishes jump-wise at a certain higher optical fre- 
quency tUop, i.e., for uj' > tOop, and the PC medium is approximately treated as a free space 



with £:(r) = 1. We choose c^op in such a way that our results are verified to be insensitive 
to perturbations. lOo^ajl-nc ~ 3.5 is chosen in our calculations. Furthermore, we distinguish 
two different types of integrals for /5(r, cu — cJj): the principal integral, when the integrand 
in Eq. (0) has a singularity, and the normal integral, otherwise. With this in hand, we 
find that the terms for j < / and for j > / in the right hand side of Eq. (jHI) contribute the 
principal and normal integrals near lOi , respectively. In order to show this clearly, we assume 
that cij = (5 + tu; is a solution of Eq. (0), and \b\ ^ cj^+i — uji, where wj+i is the closest to 
and higher than the frequency of the level /. For j < I, the integrand has a singularity due 
to 6 + ui — Uj > 0. But for j > I, the integrand has no singularity due to 6 + ui — Uj < 0. 

In PCs, the LSRF g{r,uj') displays dramatic fiuctuations when the frequency u' varies 
for a given position. As an example, we demonstrate this in Fig. 1 for an 3D inverse-opal 
PC J28l] without stacking faults 29|. Thus, the principal integral (3{r,uj) {u > 0) should 
be very sensitive to the value of uj, and the contribution to the integral comes mainly from 
the region near the frequency u. Figure 2 shows that (3{r, u) is an oscillatory function of u. 
However, for the normal integral /3(r, —uj) {uj > 0), the fiuctuation in g{r,uj') are smoothed 
out after integration, and (3{r, —u) is a slowly varying function of u, similar to the case 
of vacuum. In Fig. 3, we find a confirmation of this behavior of the function /5(r, — cj). 
Furthermore, it can be seen that in a PC the function (3{r, —u) tends to the limit value of 
that in vacuum as the frequency uj grows. Therefore, the terms with j > / in the right-hand 
side of Eq. (jH)) can be treated similar to the case of vacuum. If we consider u — ui ^ 1, 
then the PCs do not bring about appreciable changes in those terms with j > I compared 
to the case of vacuum. Therefore, Eq. (jHl) can be approximated as follows 

' ^ 27r J, {u-u,-uj')u' ^> 

Equation Q shows that, compared to the case of vacuum, inhomogeneous EM systems lead 
to an additional contribution to the Lamb shift that comes mainly from the real photon 
processes, rather than the virtual photon precesses, in contrast to the case of vacuum. 

We now apply our result to study the Lamb shift for a hydrogen atom in the inverse- 
opal PC. First, we obtain an interesting result that the PCs environment has no effect on 
the 2s state due to a2sis = 0; this result coincides with the prediction obtained earlier 
from the isotropic dispersion model [l8|. However, for the 2p state, we have Agp = and 
a2pis — 4 X 10^^. Numerical results for the Lamb shift of the 2p state are presented in Fig. 



4. We find no level splitting, which differs our finding from the prediction of the isotropic 
model |l8| . In addition, the Lamb shift depends strongly on not only the transition frequency 
but also on the atomic space position, different from dispersion models |l8|, |l9|, |20|. The 
similar properties can also be found for the 3s, 3p, and 3d states. 

Analyzing the results presented in Fig. 4, we notice that the Lamb shift can take very 
large positive or negative values and, therefore, it can be termed as a giant Lamb shift. 
Comparing the results for the PC with those for vacuum, we find that the Lamb shift may 
be enhanced in the PC by one or two orders of magnitude. Furthermore, it is significant 
to point out that the giant Lamb shift may occur for the transition frequency being either 
near or far away from the photonic band gap. The above-mentioned results are in contrast 
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2l|. In Ref. |2l|, 



to the predictions based on dramatically simplified models 
a PBG structure was simply treated as an averaged homogenous medium. This smooths 
out the contribution to the Lamb shift from real photon processes that play a key role in 
inhomogeneous systems. In the isotropic model (l 



g{r,uj) ~ (uj — ujc)^^^'^ /io, that gives an 
infinite interaction between atom and photons at the band ed ge u = Uc leading to the level 
splitting and anomalous Lamb shift. In the anisotropic model [l^, g(r,uj) ~ (cu — cUcY''^ /to, 
that leads to coupling interaction near the band edge being smaller than that in vacuum 
where q(r ,Lj) = u; it predicts much smaller Lamb shift than that in vacuum. In pseudogap 
model J20|, ^^(r, u;) ~ uj{l — hexp[{uj — cUcy/a'^]}, that gives rise to small values of the Lamb 
shift near a pseudogap. Clearly, these models lose the main physical characteristics of the 
LSRF g{r,u) in realistic PCs that may result in the giant Lamb shift and other significant 
effects. 

Based upon the position-dependent Lamb shift, we can suggest a possible experimental 
approach for verifying our theoretical predictions: if an assemble of atoms spreads randomly 
in PCs, the atoms at different positions have different values of the Lamb shift. Then the /- 
state levels of many atoms should form a /-state mini-band, similar to the velocity-dependent 
Doppler effect in atomic/molecular gases. This mini-band should be experimentally observ- 
able through the emission spectrum of these atoms. 

In conclusion, we have developed a general formalism for calculating the Lamb shift for 
multi-level atoms. It is revealed that the real photon processes play a key pole in inhomoge- 
neous dielectric structures. Our numerical results for atomic hydrogen in a 3D inverse-opal 
PC show that the Lamb shift may be enhanced remarkably by the PC environment. We 



have also predicted the existence of the Lamb shift mini-band for an assemble of atoms 
opening up possible ways for experimental observations. We believe our results provide a 
deeper insight into the theory of spontaneous emission in PCs and many applications such 
as the development of thresholdless lasers. 
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FIG. 1: Local spectral response function g(r,uj) for an atom placed at three different positions: 
ri = (0,0, 0)a, r2 = (0.34, 0,0)a, and rs = (0.24, 0.24, 0)a in the inverse-opal photonic crystal 
created by air spheres in a medium with n = 3.6 and / = 0.74; a is the lattice spacing. 
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FIG. 2: Principal integral /3(r,u;) for three different atomic positions in the photonic crystal. All 
parameters are the same as in Fig. 1. The dash line corresponds to the case of vacuum. 
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FIG. 3: Normal integral /3(r, —lv) for three different positions in the photonic crystal. All param- 
eters are the same as in Fig. 1. The dash line corresponds to the case of vacuum. 
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FIG. 4: The Lamb shift of 2s and 2p states as a function of the lattice constant a for atomic 
hydrogen located at three different positions in the photonic crystal. The Lamb shift of 2s state 
(the dashed lines) is the same as the value in vacuum. All parameters are the same as in Fig. 1. 
The frequency region between two dashed-dot lines is the photonic band gap. 
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